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Abstract 

We analyze sum rules for the T system with resummation of threshold effects 
on the basis of the nonrelativistic Coulomb approximation. We find for the pole 
mass of the bottom quark mj, = 4.75 ± 0.04 GeV and for the strong coupling 
constant a s (Mz) = 0.118 ± 0.006. The origin of the contradiction between two 
recent estimates obtained within the same formal framework is clarified. 



1 Introduction 



The knowledge of precise numerical values of the coupling constant and quark masses 
is important for obtaining accurate predictions within QCD. The technique of sum 
rules is extensively used both for the determination of these parameters from exper- 
imental data on low lying resonances and for the calculation of properties of higher 
excited states [l], 0. The T system is a rather clean place to obtain the values of 
the strong coupling constant a s and also of the 6-quark pole mass m& which is a key 
parameter for the rich B-meson physics. In sum rules for the T system there are 
some subtleties that make them interesting by themselves, the main one being the 
importance of Coulomb effects and the applicability of the nonrelativistic approxi- 
mation as a rather good starting point for the expansion in the coupling constant 
a s and the 6-quark velocity v which are both small parameters of the problem ||. 
The expressions for the first few moments of the spectral density are now available 
in ordinary perturbation theory (PT) with a 2 s accuracy [|J, however, they cannot 
be used in theoretical formulas for sum rules directly because the spectrum is well 
known experimentally only for energies close to threshold due to existence of sharp 
resonances while the contribution of the continuum to these low moments is large in 
comparison with the resonance contribution. For higher moments that are saturated 
by these low energy resonances with better accuracy the ordinary PT expansion in the 
coupling constant is not applicable (the coefficients of the expansion grow fast with 
the order of the moment) and a resummation is required. In this paper we analyze 
the sum rules for the T system with resummation of Coulomb effects that are largely 
responsible for the growth of PT coefficients to extract the values of a s and m&. We 
consistently take into account the first order corrections in a s and give an estimate 
of the magnitude of the high order corrections. 

The analysis is based on the same formal approach which has already been applied 
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in refs. pi || along different lines and with contradictory results. In ref. || the next- 
to-leading approximation in a s both to the high energy part of the correlator and 
to the Coulomb potential was used while in ref. || further terms of the perturbative 
expansion in a s have been partly included as well. This prevents the direct comparison 
of the results of those two papers though both claim that the corrections are small and 
under control. However, the numerical results do not overlap within the respective 
error bars and the clarification of the origin of this difference is in order. We find 
that the difference is to large extent due to neglecting in ref. |5| the contribution of 
Coulomb poles to the moments. We demonstrate that after including the contribution 
of the Coulomb poles the results of two papers can be made completely compatible 
by the corresponding choice of the scale in the nonrelativistic Coulomb potential - 
the low one (// ~ 1 GeV) for the case of ref. || and of order of the 6-quark mass for 
ref. §. 

The result of ref. || for the value of a s [Mz) is furthermore in contradiction with 
the LEP data 0. The main uncertainty enters the theoretical predictions through 
the higher order a s corrections which become uncontrollable at the low normalization 
point used in ref. ||. The uncertainties of the predictions given in ref. || are there- 
fore strongly underestimated explaning thus the discrepancy with the LEP value for 
a s (M z ). 

The paper is organized as follows. The stage is set in Sect. 2 where our notations 
are fixed. The numerical analysis of sum rules is done in Sect. 3 where the procedure 
of optimization for the nonrelativistic Coulomb Green function is described. In Sect. 4 
we discuss the difference between the present analysis and those of refs. ||, ||. Sect. 5 
contains our conclusions. Some explicit formulas are collected in the Appendix. 
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2 Sum rules and nonrelativistic approximation 

We perform our analysis by computing the moments directly instead of using a repre- 
sentation for them through the propagation function of the nonrelativistic approxima- 
tion in Euclidean time as done in ref . || . The difference between these two approaches 
is briefly discussed later. 

The moments of the spectral density are defined to be dimensionless and are given 
by the following expression 

R(s)ds 



1 2tt 2 d n 



(Amir 



CO 



S n+1 



where 

R(s) = 127rlmn(s + ie) 
and n(s) is the polarization function of the 6-quark vector current 

(qtfu ~ 9„q 2 ) I% 2 ) = % J dxe^(0\TUx)U0)\0), 

jfi — b'y^b. The assumption of quark-hadron duality relates the imaginary part of 
the polarization function R(s) to the experimentally measurable quantity Rb(s) = 
a(e + e~ — > hadrons^) /cr(e + e~ — > n + n~) by the sum rules 

r°° R(s)ds (Am 2 b ) n f°° R b ( s )ds 



re 

M n = (AmlT / 
Jo 



s n+1 Ql Jo s n+1 

with Qb = —1/3. 

The experimental values of the moments are obtained by using information on 

the first six T resonances. Their leptonic widths and masses (k — 1 ... 6) are 

known rather precisely |7j 

(Amir f ^ ^ T fc too R b (s)\ 
n ~ Ql W QED (m b )hM! n+1+ Js dS s^) W 

while the last term is a continuum contribution. Beyond the resonance region, for 

energies larger than s ~ (11.2 GeV) 2 , the spectral density Rb(s) can be approximated 



by the ordinary PT expression for the theoretical spectral density R(s). The influence 
of the continuum on high moments is almost negligible numerically and in any case 
under strict control. Here an ED (mb) = 1.07a 2 @. The theoretical expression of 
the polarization function in the nonrelativistic approximation that is justified for the 
energy region near the threshold is given by || 



where Cp = 4/3. The energy of the nonrelativistic Coulomb bb system counted from 
the threshold is E = —k 2 /mf,, and there are several ways of introducing k through 
the relativistic parameter of energy square s 



which all are equivalent in nonrelativistic approximation near the threshold. For 
E < (s < Ami) the parameter k is real while for E > (s > Ami) the parameter k is 
purely imaginary. The quantity G(0, 0, E) is the nonrelativistic Green function (GF) 
of the bb system that can be found exactly for Coulomb potential which is supposed 
to dominate the whole QCD interaction in the energy region of interest. The ct s (n) 
correction in eq. (|2|) (so called perturbative factor N pt ) reflects the contribution of 
high energies to the polarization function and has as natural scale /i ~ m^. We use 
the MS scheme of renormalization throughout the paper and m^ denotes the pole 
mass of the b quark, as already mentioned in the Introduction. The PT expression 
for the nonrelativistic potential is known up to the second order in a s 



where Vc is the Coulomb potential, Vq = —Cpa s /r and explicit expressions for V\,t, 








are given in the Appendix (eq. (0)). The corrections to the Coulomb potential in 
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eq. (|]) generate corrections to the Coulomb GF Gc- The full nonrelativistic GF in 
condensed notation (see Appendix for more details) is given by q 

G = G c + AG = G c - G c ® AV ® G c = G c + G 1 + G 2 . (5) 

where Gc is a Coulomb GF and G\ and G2 correspond to V\ and V2 contributions 
(see eqs. (g]| |0|, 0) of Appendix). 



In principle one can develop a regular expansion for the polarization operator near 
the threshold using Bethe-Salpeter equation for four-point amplitude using Barbieri- 
Remiddi exact solution |10| as a starting point. This is however a rather difficult 
technical problem. Meanwhile up to the first order in a s eqs. (Q, ||) give a complete 
result. The self-consistent calculation of the expansion of the polarization function 
beyond the first order is not available, so we will use the G2 term of eq. fl5|) not for the 
regular computation but rather for estimating the typical high order a s contributions 
to the sum rules. 



3 Numerical analysis of sum rules 

We consider two methods to compute the nonrelativistic GF when higher order per- 
turbative corrections to the potential are included. The first one is the direct computa- 
tion when the splitting into the leading term (Coulomb) and corrections is determined 
by the ordering in the MS scheme, i.e. we count powers of a s . This method in fact 
fails because the expansion turns out to be badly divergent. The second methods is 
the optimized computation when the splitting of the whole potential into the leading 
term (Coulomb) and corrections is determined by the physics of the problem rather 
than the artificial MS charge. This technique gives solid results when all expansions 
(for the effective charge and for the nonrelativistic GF) converge well for the moments 
that are computed. 
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3.1 Direct computation 

To compute the theoretical part of sum rules we use the approximation (|2]) for the 
polarization function. As a function of the nonrelativistic variable E the exact non- 
relativistic Coulomb Green function of the bb system G(0, 0, E) has its poles below 
threshold for real E < (s < Ami) corresponding to the bound states, and a cut 
above the threshold at real E > (s > 4ml). Accordingly, the theoretical moments 
are split into 

M n = V n + C n 

where V n is the pole contribution from the energy region below the threshold 

r 4 

V n = (Amir / 
Jo 

and C n is the continuum contribution from the energy region above the threshold 

R(s)ds 
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v c 

' n 


V 1 

' n 


V 2 

' n 


10 


0.16017 


0.23375 


0.25756 


12 


0.16631 


0.25180 


0.27743 


14 


0.17271 


0.27101 


0.29859 


16 


0.17939 


0.29145 


0.32111 


18 


0.18635 


0.31319 


0.34507 


20 


0.19360 


0.33630 


0.37057 



Table 1: The 0th, 1st and 2nd order pole contributions to the moments 
Then the total moments are represented by the sum 

M n = M° + Mi + M 2 n (6) 
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where is the contribution of a pure Coulomb part, A4 l n are the contributions 
due to Gi terms in eq. (|). The same representation is used for V n and C n parts 
separately and also for H(s) and R(s), for instance, V n = V% + + V% anc ^ so 
on. An example of the breakdown of the total numerical values of the moments into 
particular contributions due to Gi terms for \i = = 4.75 GeV, a s (M z ) = 0.118, 
that are chosen as our typical values (see below), is given in Tables [p. |2|. 



n 




C 1 

Si 




10 


0.26388 


-0.02508 


0.00320 


12 


0.22195 


-0.02055 


0.00252 


14 


0.19059 


-0.01724 


0.00200 


16 


0.16652 


-0.01477 


0.00159 


18 


0.14763 


-0.01288 


0.00125 


20 


0.13247 


-0.01139 


0.00099 



Table 2: The 0th, 1st and 2nd order continuum contributions to the moments 

From the Tables 1,2 one sees that the expansion (|5p does not work for the pole 
contribution even at relatively large scale (/i = m^), i.e. next terms are larger than the 
lower order terms. In Voloshin's paper these corrections were treated by a redefinition 
of the scale in the leading order approximation. This resulted in the introduction of 
the function ^(7) that determines the appropriate scale for a given moment. Because 
only the first order corrections (Vi from Eq. (||)) were taken into account (the contri- 
butions due to V2 were not considered) such a redefinition of the scale in the leading 
order approximation does allowed to absorb them completely into the Coulomb Green 
function requiring G\ = 0. However, this leads to an extremely low normalization 
point where the ordinary PT for the potential breaks down and high order corrections 
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become uncontrollable (see Sect. 4 for the details). 
3.2 Optimized computation 

Now we turn to our method of computating the nonrelativistic Green function with 
optimization of contributions due to higher order corrections to the potential. In con- 
trast to QED where the Coulomb term exhausts the potential in the nonrelativistic 
approximation and hence generates the full Green function, in QCD the potential 
itself is a PT expansion in a s . Due to the presence of log(fir) terms in this expan- 
sion the nonrelativistic Green function cannot be obtained in closed form and some 
approximations have to be employed. The perturbative expansion in the corrections 
to the potential looks natural for this purpose. Now, for computation of the best 
approximation for the nonrelativistic GF (|5]) we apply the technique of optimization 
of the expansion @ with regard to higher order corrections using the freedom of 
redefinition of the leading order approximation by the proper choice of the Coulomb 
coupling constant ||11|| . It is clear that a s in the MS scheme is not necessarily the best 
choice for the coupling if the Coulomb solution is used as a leading approximation for 
the nonrelativistic GF. To improve the situation we introduce an effective Coulomb 
charge by adding a part of the higher order constant contribution to the leading term. 
We parametrize this part with a parameter t. To some extent this is reminiscent of the 
simple change of the scale /i, though in reality it is a more sophisticated procedure. 
Its intension is to minimize (or optimize) the contribution of the logarithmic terms in 
the potential that cannot be treated exactly. The fi dependence is still present and 
controlled by the ordinary renormalization group invariance. Moreover, in contrast to 
ref. H where different normalization points are used for the analysis of the soft and 
hard corrections, in our approach we may use a unique expansion parameter a s (m(,) 
and the corrections related to the matching are absent. Thus, the Coulomb part of 
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the GF is determined by an effective Coulomb parameter otc(t) 

V c = - Cpac ^ (7) 
r 

with acif) being a function of a s 

a c (t)=a s (l + ^(l-t)C^. (8) 

The explicit formula for such a decomposition to one loop (first order approximation 
for the potential) reads 

C F a s ( ^ , a s x x 



V(r) = -^(^l + ^^ + ^lnM) 

^ f l + ?(1 - t)Ci) - + C\ In 



r \ 47r / r An 

= Vc+^-V, (9) 

47T 

with C\ given in the Appendix. Now 

471 47T 

where the effective Coulomb term Vc is given by (|7]) with the effective charge (§) 
and Vi is considered as a correction. Thus we have added a part of the constant 
of the second order term (Vi) to the Coulomb charge of the leading approximation. 
Developing the expansion around the Coulomb GF corresponding to Vc we minimize 
corrections by choosing the proper numerical value of the parameter t. At the final 
stage we have to check that both expansions (for the nonrelativistic GF around the 
Coulomb approximation and for the effective charge around the a s ) converge well. 
This happens to be the case for the considered system. Otherwise one would have 
to conclude that the problem cannot be treated perturbatively, i.e. one cannot use 
the PT expression for the potential to generate the nonrelativistic approximation. 
Technically the method results in substituting the combination Cfol s by Cpotc in Gc 
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and G\. Then we choose t(n) (or ao) such that M. l n = to get the best convergence 
of the series (0). 

In the analysis we use moments with 10 < n < 20. For this range of moments 
the contribution of the continuum (not well known experimentally) is sufficiently 
suppressed for the results being practically independent of s and, at the same time, 
the nonperturbative power corrections due to the gluonic condensate are suppressed 
too||. We use the normalization point \i = m b and take into account the first order 
a s corrections to the polarization function (0). 

The result of the fit a s (m b ) = 0.22 ± 0.02 or a s (M z ) = 0.118 ± 0.006 is in a 
good agreement with other estimates 0. The sum rules are much more sensitive 
to the 6-quark mass than to the strong coupling constant so it is instructive to fix 
a s (Mz) = 0.118 to the "world average" value and then to extract rrib{n). The 
results for rrib(n) are given in Table [|, together with the individual moments and the 
respective values of the effective charge. 



n 


m b (GeV) 




c 


«c/«s 


10 


4.730 


0.4143 


0.3308 


1.306 


12 


4.742 


0.4723 


0.2827 


1.328 


14 


4.752 


0.5391 


0.2481 


1.349 


16 


4.760 


0.6028 


0.2174 


1.362 


18 


4.767 


0.6786 


0.1949 


1.376 


20 


4.774 


0.7642 


0.1780 


1.389 



Table 3: The n distribution of m b , V n , C n and ac/ot s for \i = m b , a s (M z ) = 0.118 

Our final estimate of the 6-quark mass is m b = 4.75 GeV, the average over the range 
n = 10 — 20. As we see, the pole contribution exceeds the one from the continuum and 
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obviously cannot be neglected as it was done in ref . [0 . What is the uncertainty of this 
result? The total error from the specific choice of k(s) in (|3]) and from the numerical 
value for so is about ±0.1%. However, the main uncertainty originates from higher 
orders of PT for the potential @ or from the expansion of the GF in eq. (|5]). This 
can be estimated from the n distribution (Table |3|) or from the /i dependence of the 
results which formally must be canceled by the higher order contributions (Table 



li(GeV) 


m 6 (GeV) 


N pt 


a c /a s 




m b 


4.752 


0.631 


1.349 


0.218 


4.0 


4.763 


0.608 


1.254 


0.231 


3.0 


4.782 


0.564 


1.245 


0.257 


2.0 


4.799 


0.478 


1.118 


0.308 


1.2 


4.767 


0.273 


0.879 


0.428 



Table 4: The fi dependence of nib, N pt , otcl&s and a s for n = 14, a s {Mz) = 0.118 

Note that at the scale // ~ nib both the hard gluon corrections which are parametrized 
by the factor 

N pt = l-AC F ^ (10) 

7T 

and the soft gluon corrections parametrized by the ratio a.c/ot s are of the same order 
~ 35% so this seems to be an optimal normalization point. For lower scale the hard 
corrections become large while for higher scale the same is true for the soft corrections. 
These latter corrections can be also estimated directly by adding the G2 part of the 
nonrelativistic GF. The relative weight of a s corrections can be found from Table [5[ 
In the Coulomb approximation the hard gluon corrections, i. e. the perturbative factor 
in eq. (Q) is neglected while the 1st order corrections include this factor and the G± 
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part of the nonrelativistic GF. In the second order estimate we take into account the 
G 2 term but do not use the optimization procedure. Thus collecting all the data we 
estimate the high order corrections to the 6-quark mass to be about 1%. Our final 
value is m b = 4.75 ± 0.04 GeV. 



Order of approximation 


m b (GeV) 


pure Coulomb approximation 

total 1st order a s corrections included 

2nd order a s corrections to the potential included 


4.700 
4.752 
4.776 



Table 5: The convergence of perturbation theory for m& at n = 14, /i = m^, a s [Mz) — 
0.118 

Let us emphasize that we have used the a 2 s corrections to the potential only 
as an estimate of the typical size of all a 2 s corrections. The final answer in the next- 
to-leading order can be obtained only after all a 2 s corrections (including relativistic 
ones) to the polarization function in the threshold region are available. 

4 Comparison with the previous analyses 

In ref. [||] the low normalization point fi = 1 GeV is used which is effectively a way of 
reshuffling the PT series in higher orders. The hard gluon correction was separated 
and taken at fi = 0.632m;,. Using these prescriptions we reproduce the results of ref. 
|m i = 4.827 ± 0.007 GeV, a s (M z ) = 0.109 ± 0.001 within our approach. 

The precision of these results is however strongly overestimated: no variation due 
to higher order corrections (which at this level reduces effectively to the change of 
the parameter /i) is allowed. The quoted errors just reflect the difference between 
different moments while the /i dependence is fixed from the very beginning. Also, 
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although the correction induced by using different normalization points for hard and 
soft regions respectively is formally of the next (second) order in a s , they can be very 
essential because of the low normalization point. This can be seen clearly from the 
fact that at fi ~ 1 GeV the first order hard gluon corrections are about 80% so that 
PT does not converge. On the other hand we found that the second order correction 
to the moments due to Vi term of the potential (§]) exceeds 100% at /i ~ 1 GeV so 
the first order PT can be hardly used for the analysis of the soft gluon corrections at 
this scale though the first order corrections are minimized (equal to zero). The error 
bars given in ref. || do not account for high order corrections. This is the reason for 
the discrepancy between the result of ref. || and LEP data on a s [Mz) 0- 

Speaking of uncertainties we would like to touch the point of using the moments 
in relativistic representation against the Laplace transform of the nonrelativistic GF. 
We mentioned this problem in the beginning of Sect. 2 when explaining details of 
our treatment of nonrelativistic approximation. The parametric estimate of 0(l/n) 
for difference between moments and Laplas transform in ref. || is not correct for an 
arbitrary large n. Let us demonstrate it in more detail now. Two formulas for the 
moments 



that are used and where s = (2m + E) 2 are not equivalent up to 1/n terms for large 
n for any function R(s) but only for rather smooth ones. In fact, the dominant 
contribution to the theoretical value of the moments is due to Coulomb bound states 
which lead to ^-function singularities in R(s). So, one has to check the equality 



for R(s) = Aml5(s — ttiq) with mc = 2mf>(l + A) where A ~ a 2 , A < 0. Then one 
has 




vs 





(11) 
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while the right hand side of eq. (|TT| ) is 

dE e- 2An 



f°° -JL n dti e 
/ e m 6 R(E) — = - 
Jo rrih 1 



m b 1 + A 
or 

(TT^ vs (12) 

and these two quantities are not related by the estimate flTID, i.e. with an accuracy 
of 0(l/n) uniformly in n and A. Indeed for large n and fixed negative A one has 



e 2Ar, 



;i + a) 



2n+l 



0*l + O[~). 



n 



n— >oo 

They are close to each other with the stated accuracy only for \A\ 2 n << 1 that can 
be easily checked by expanding both quantities in a Taylor series for small |A|. For 
a smooth function, however, the estimate (|TT| ) is valid. For instance, for R(s) = 
(s/Am 2 b )6(s - Am 2 b ) = + E/2m b ) 2 9(E) one has 

9x „ f°° R(s)ds , 9N „ i r°° sds 111 
4^6 / ' = 4m "- 1 / — — - = = - + — + ... 13 

V b> Jiml S n+1 K U JAmi S n+1 n - 1 U U 2 V ' 



b 



while 



r°° -B_ nn ,„ s dE r°° -JL n ( E E 2 \ dE 1 1 1 
jo nib Jo \ vrib 4m^ J mb n nr In 6 



(14) 



in accordance with eq. fllTf ) and eq. (TTT|) . This observation in fact gives an additional 
constraint from above on possible n. Meanwhile in ref. the upper bound on allowed 
values of n was connected only with the contribution of the gluon condensate. We, 
however, should stress that in ref. || the region of the parameters with |A|n ~ 1 
(a 2 n ~ 1) was considered where the difference between two expressions in eq. ( |12|) 
can be estimated numerically as 0(l/n) or 0(a 2 ). 

The results of ref. m b = 4.60±0.02 GeV, a s {M z ) = 0.119±0.008 are reproduced 
as soon as the pole term is omitted. This is the reason for having smaller value of 
m b . In fact, omitting the pole contribution leads to a huge (about 50%) change of the 
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moments (see Table |^) but only to ~ 3% variation of the 6-quark mass. Then the way 
the potential is treated in ref. B seems not to be correct to us. A simple use of the 



running coupling constant ay{q) with q = y s — Ami which is defined in all orders of 
PT as V(r) = —Cpay(r)/r instead of a s in the leading order nonrelativistic result 
(Sommerfeld factor) || is not consistent with the expansion of the GF (|5|) that is 
generated by PT of the potential (§) though it does catch a part of the contribution. 
This can be directly found by comparing the data listed in Table ||] and Table |6| which 
contains the continuum contributions Z n to the moments computed according to the 
prescription of ref. |J by inserting the running coupling ay{q) into the Sommerfeld 
factor. Note that in our approach this procedure corresponds to inserting ay{q) 
instead of a s to the Coulomb GF Gc- In Table [| the quantities Z % n are defined as 
follows 

/•OO A e 

Zi = (4m 2 b ) n —tR C (s)\ Wm (15) 

J Am 2 S n+1 n as—a y v '(q) v ' 

b 

where ay is a series for ay in a s up to the zth order 

4 0) = a s , c® (q) = a s (l + ^ L - 2/3 In ( £\\ V ... 



( a v (?) an d numerical values of parameters can be directly found from eqs. (|16| , 
in the Appendix). The quantity R c (s) for s > Ami reads 

, 1 8tt 9 AT / 7rCVa s \ 



ml ' 'l ' \l — exp(—7rCFa s /v) J ' 

■u being the 6-quark velocity. Note that in eq. (15) a s is not substituted by ay in the 



perturbative factor N pt (PlJfUl). Thus the quantities and Z\ should be compared 
with C% + C\ and C% + C\ + C% respectively. From Tables ||, || one finds that within 
the approach of ref. |J the corrections to the leading order continuum contribution 
are about 2-4 times underestimated. In any case, the main difference between the 
present analysis and that of ref. H is just the neglect of the pole contribution. 
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n 




Z 1 


Z 2 


10 


0.26388 


0.25534 


0.25273 


12 


0.22195 


0.21564 


0.21380 


14 


0.19059 


0.18582 


0.18450 


16 


0.16652 


0.16285 


0.16188 


18 


0.14763 


0.14475 


0.14403 


20 


0.13247 


0.13018 


0.12964 



Table 6: The corrections to the continuum contribution computed by inserting the 
running coupling ay[k 2 ) into the Sommerfeld factor for /i = m&, oc s (Mz) = 0.118 

Furthermore the results in ref. || exhibit a much weaker /i dependence and the 
contribution of high orders of PT is claimed to be much smaller. The reason is that 
the main corrections and the strong /i dependence originate from the pole contribution 
which is omitted in ref. 0. 

5 Conclusion. 

We reanalyzed the sum rules for T system and obtained new numerical estimates 
of the strong coupling constant and the bottom quark pole mass. Our results are 
different from those two existing papers on this subject and lie outside error bars 
given by the authors. Our results and the reasons for the discrepancy are: 

• We found m b = 4.75 ± 0.04 GeV and a s (M z ) = 0.118 ± 0.006 consistent with 
LEP data. The results are insensitive to the choice of k in eq. (|||), to s , and 
whether resumed or plain PT is used for integration beyond Sq. The main 
theoretical uncertainty of these values comes from the higher order corrections 
which induce an n and /i dependence of the result. It is explicitly shown that 
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the part of the corrections connected to the potential is large and leads to less 
accurate estimates of the mass and the coupling constant than claimed before. 
Thus the main problem at the present stage of calculation is the consistent 
treatment of the next-to-leading a s corrections. 

• The pole contribution to the theoretical spectral density is omitted in ref. ||. 
This results in the a significant difference in the determination of mj, (which we 
consider as erroneous) and leads to the impressive stability against a variation 
of \i (or against inclusion of higher orders of PT) because it is the pole part 
that is strongly \i dependent. 

• In the paper [E] the low normalization scale is used where the validity of the 
expansion around Coulomb solution due to huge high orders a s corrections 
is questionable. The uncertainty of the results is strongly underestimated in 
ref. §. 
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Appendix 

The corrections to the Coulomb potential read || 

V c {r) = , 

V 1 (r) = V c (r)(C 1 +C 1 l \n(r f ,)), 
V 2 (r) = V c (r)(C 2 + C?ln(r/i) + Cfln 2 ^)), 

where 

Co = ai + Zfio^E, 
19 



V 



Ci =[y + A ^J Po + 2(01 + WoaihE + a 2 , 
C\ = 2ft, 
C\ = 2(& + 2/3 a 1 ) + 



C£ = 4$, (17) 



-(y-16C(3))^T F n / +( Y T F n / )' 



4 1 
Ca = 3, Op — -, Tp = -, 

7e = 0.57. . . is the Euler constant, ((x) is the Riemann (-function, and rif is the 
number of light flavors. 

The nonrelativistic Coulomb GF reads 131 



47T \ Cpa s m b \Cpa s m b 1 ^ 



m) 



m=0 



where 



F(m) = L . (19) 

. n C F a s m b m + 1 

m + 1 : 

2fc 

Evaluation of the corrections to the Coulomb GF is straitforward but rather compli- 
cated. In the calculation we used the representation of the Coulomb GF Gc(x,y, k) 
as a sum of Laguerre polynomials |J . The result reads 

0,(0,0, k) = (£ F\m){m + 1) (c\ - In ^ C? + 2)C^ - 

oo m-1 _,_ i 

2 E E 

m =l n =0 " 
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2 E F ( m ) ((Co " ^ (j) C^J - {2 lE + Mm+ 1))C7^ - 
G2 (0,0,^( S )^(| / V)(( m+1 )(^-ln(f)^ + 



In 2 Cjj + (m + l)*,(m + 2) fa 2 - 2 In ^ C 2 j + /(m)C|j + 
WW ("^ (<S - 2 ln (f) Cl) + J(n ,,n)C 2 2 ) + 

2| /W ((^-lng) C ? + .n 2 (f) C2 2 )- 
(2 7 e + *!(m + 1)) fa 2 - 2 In (jj) C 2 ") + A'(m)C 2 2 j - 



where 



/ 7T 2 2 

I(m) = ( m + 1) # 2 (m + 2) - ^ 2 (m + 2) + — 



3 (m + l) 2 / 

2(tt 1 (m + l)+ 7j5 ), 

n + 1 / 1 \ 

J(m, n) = 2 *i(m - ra) + 2-f E ) + 

m — n \ n + 1 / 

777- H - 1 

2 (*i(m - n + 1) - *i(m + 1)), 

m — n 

KT(m) = 2(*i(m + 1) + 1e? + ^2(m + 1) - ^ 2 (m + 1) + 2 7 |, 
and r(x) is the Euler T-function. 
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